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Moving point charge as a soliton in nonlinear electrodynamics 
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The field of a moving pointlike charge is determined in nonlinear local electrodynamics. As a model 
Lagrangian for the latter we take the one whose nonlinearity is the Euler-Heisenberg Lagrangian of 
quantum electrodynamics truncated at the leading term of its expansion in powers of the first field 
invariant. The total energy of the field produced by a point charge is finite in that model; thereby its 
field configuration is a soliton. We define a finite energy-momentum vector of this field configuration 
to demonstrate that its components satisfy the standard mechanical relation characteristic of a free 
moving massive particle. 
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I. INTRODUCTION 

In this paper we are studying the field of a uniformly moving charge in nonlinear electro¬ 
dynamics, in other words we are dealing with nonlinear extension of the Lienard-Wiechert 
potentials. 

The problem of the field caused by a moving charge, besides belonging to fundamental 
problems of electrodynamics, is also of certain practical importance as applied to charged 
particle beams in accelerators [l[. Its nonlinear extension may have an effect in small- 
impact-parameter scattering at high energies where close vicinities of the charge come into 

play- 

We take the simplest version of nonlinear electrodynamics , namely the one that results 
from keeping only the leading term in the expansion of the Heisenberg-Euler action Q of 
quantum electrodynamics (QED) in powers of the field invariant 5'(cc) — 

A^), with the second invariant 0(a;) =\Ffj^^F^'' = B • E set equal to zero. (Remind that, the 
Heisenberg-Euler action reflects the nonlinearity due to interaction between electromagnetic 
fields stemming from their quantum nature, and it is local in the sense that dependence on 
space-time derivatives of the fields present in the full effective action of QED is disregarded 
from it.) This nonlinear model may be thought of as an extension of the established theory 
of electromagnetism to the close vicinity of a point charge compatible with what that theory 
describes at larger distances. When treated seriously, the above truncated model results in 
the electric field behavior near the point charge less singular than the Coulomb law, which 
provides convergence of the field-energy integral 0 . This means that the field configuration 
produced by a point charge may be referred to as a soliton. In Section IIIII we present the 
electric and magnetic fields of this soliton when it moves with a constant speed. In Section HVl 
we define the energy-momentum of this moving soliton, which is just the energy-momentum 
vector of a moving particle with its mass defined by an integral of its field. It is proportional 
to the charge squared. 

The following remark is in order. Our approach to considering the nonlinear problem 
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includes as its starting point a covariant representation for the current, the field, the energy- 
momentum tensor and vector. In the trivial, linear limit, this approach results in totally 
reproducing the field of a charge moving with a constant speed as given in the Landau 
and Lifshitz text-book [^. This seemingly undistinguished fact is worth stressing, however, 
because, strange though it may seem, one faces a discord in the present-day literature, many 
basic text-books included, concerning this fundamental point (see the analysis in Q). The 
dispute is fueled by the expression for the field of a moving charge given without a detailed 
derivation in Feynman’s lectures @ , that may seem to be different from that of Ref. Q , the 
more so that apparently different plwsical ideas are layed into it. This fact even called into 
being an experimental verification [l[ that resulted in favor of [J . For this reason we found 
it not out of place to analyse the result of 0in a separate forthcoming publication Q to 
reveal its complete coincidence with that of [J both for accelerated motion of the charge 
and its motion with a constant speed. 


II. MINIMALLY NONLINEAR MODEL 


The Lagrangian of a minimally nonlinear electrodynamics is defined as 

L(x) = -^(x) + |(J(x)f, (1) 

where S(x) makes up (with the reversed sign) the Lagrangian density of the standard linear 
Maxwell electrodynamics, while the second term in © is the quartic in the field-strength 
addition to it. The field-strength tensor is related to the four-vector potential A^(x) as 
= d^A‘'(x) — d''A^(x), where Hence, the first pair of the Maxwell equations, 

[V X B] = 0 and [V x E] -|-9°B = 0, with the electric and magnetic field strengths Ei = 
and Bi = remains standard. The self-coupling constant cr is presumably small 

enough. It has the dimensionality of inverse fourth power of mass [l^ . The four-dimensional 
scalar product is {ux) = UqXo — {u ■ x) , x^ = Xq — |xf. We use Lorentz-Heaviside system 
of units. If the quartic correction in o is understood as the lowest term in the expansion 
of the Heisenberg-Euler Lagrangian in powers of S’, with 0 = 0, the value 0 of the photon 
selfcoupling constant is cr = e^/i/(457r^m'^c^), where e and m are electron charge and mass, 
respectively. 

The action includes the selfinteraction of the electromagnetic field and interaction with a 
current j'^, produced by a moving charge. It is 

5[A] = -i Ji-^A^r - Lix))d^x. 

The second pair of Maxwell equations should be calculated in assumption that the trajectory 
of the charge is fixed, and that its variation, and, thereby, the variation of the current, are 
zero. Thus, the variation of the action takes the form 


(55' = -- 


-J^SA^ 

c 


{l-a^{x))6d{x) 


(in. 


After integrating by parts and taking into account that there is no field at infinities the 
latter equation takes the form 


6S = -- 




SAfidFl 


( 3 ) 


As always, variations should turn to zero for any (5A^, and we come to the nonlinear field 
equation 


d 


[{l-a^{x))Fn = 


( 4 ) 
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III. FIELDS OF A MOVING CHARGE 

A. Current 

In this paper we restrict ourselves only to those currents that correspond to a pointlike 
or any other charge, which is at rest in the origin in a certain Lorentz reference frame (the 
rest frame) and is distributed in a spherically symmetric and time-independent way in that 
frame. The charge moves as a whole with the speed — v, z; < 1, along axis 1, Vi = Suv, in 
the inertial frame that moves relative to the rest frame with the 4-velocity which is the 

1/2 

unit four-vector, v? = l,with the components mq = Ui = 7 ^, where 7 = (l — 

The four-current is 


= u^p{x), (5) 

where p{x) is the Lorentz scalar defined as the time-independent charge density in its rest 
frame. This implies that under Lorentz boosts and spacial rotations only the argument of 
p is transformed. 

To make sure of the validity of the four-vector representation ([5]) let us first define the 
current in the rest frame as = cp{x.'), = 0. All quantities relating to that frame 

are marked with primes throughout. Applying the Lorentz transformation to this current 

jo = 7 (jo + ^j'l) > ji = 7 [j'l + ^j'o) , 32 = 42, 43 = 43 , (6) 

with the account of the vanishing of its spacial component in the rest frame j'l = 0, and 
applying the inverse transformation to the coordinates 

2^0 = 7 (a^o - = 7 , X 2 =X 2 , x'^ = X 3 . (7) 

we get for the moving charge m j°{x) = 7 cp((xi —^xo) 7 , a; 2 , ^ 3 ), j(x) = 

V 7 p((xi —-^xq) 7 ,X 2 ,X 3 ). This agrees with Eq. ([S]). For the pointlike charge e we have 

j'°{x') = cp(x') = ce<5^(x'), j'(x') = 0 

4 °(x) = 7ce(5((xi--xo ) 7)i5(x2)(5(x3) = ce 6 {xi--xo) 6 {x 2 )S{x 3 ) = 'yceS^{x'). 

\ c / c 

j(x) = -vjeS ^^xi-^xo^ 7 ^ 6 {x 2 )S{x 3 ) = veS ^xi - ^xq^ S{x 2 )S{x 3 ) = ( 8 ) 

Note that our final expressions for the current components differ by the Lorentz-factor 7 
from the corresponding expressions in Ref. 0, because the charge density p{x) is defined 
there as the zero-component of the four-current, and not as a Lorentz scalar, like here. The 
both ways are equivalent, of course. 

Later we shall see that the current ([5]) is reproduced as the right-hand side of the field 
equation (jj]) for a moving point charge in nonlinear, as well as linear, electrodynamics, and 
also in the weak continuity equation (1371) for the energy-momentum tensor. 


B. Govariant ansatz 

Once the motion of the charge is given, there are no other vectors in the problem besides 
Ufj, and the four-coordinate of the observation point Xp,, which is the argument of the dif¬ 
ferential equations ®. Therefore, the potential produced by the charge may have the only 
representation 


= u^fi{xu, x^) -I- x^f2{xu, x^). 


( 9 ) 
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where /land /2 are functions of the indicated scalars. Then the field-strength tensor resulting 
from this potential is 


with 


= {u^Xu - UyXfj)f{xU,x'^) 


f{xu,x^) 


dh _ g dfi 
d[xu) d(x^)' 


( 10 ) 

( 11 ) 


The requirement that in the rest frame the (electric) field Tpi b® independent of time Xg 
implies that / should be a function of the combination 

= (xuf - x^ (12) 


of its argument. So, by definition, the invariant W is the distance from the observation 
point to the charge in the rest frame of the latter W = |x'| = r. (In the rest frame it is 
evident that Eq. m defines a positive quantity IT = |x'| = r. As long as is a Lorentz 
invariant it remains positive in any frame.) Therefore, in what follows we shall keep to the 
representation 


f{xu,x^) = g (IT^) . 

It can be directly verified that the field tensor resulting from m and m 

— {u^Xjj UjjX^')g (yV ) 

satisfies the first pair of the Maxwell equations (the Bianchi identities) 

9F^x _ 0123 _ 1 

dxP ’ 

However, it is easier to see this if we note that by choosing /2 = 0 in and 


(13) 

(14) 

(15) 

taking /i as 


MW^) = l 19iW^)dW^ 

we determine the vector-potential generating the field ([T4|) . (IT^ 

Ap = Up/i(lT2) 

thus guaranteeing the fulfillment of ((TS|) . In the rest frame this potential is the 3-scalar A'g = 
/i(x'^). A' = 0. Setting /2 = 0 is not the only way to fix the vector potential generating 
Eq. (fT4)) . Another choice of the potential admitted within the gauge arbitrariness may be, 
for instance, 

Ap = (up [ux] - Xfi) (ux) g{W^). 

This potential satisfies the Lorentz-invariant gauge condition (Au) = 0. 

In what follows we concentrate in finding a solution to the second pair of the nonlinear 
Maxwell equations (|4]) using the form (IT4l) as an ansatz for searching for it. 


1. Linear limit 


Let us first note that the covariant extension of the Coulomb field produced by a moving 
point charge e, which is at rest in the origin in the rest frame, i.e., by the one, whose world 
line passes through the point xg = Xi = 0, is: 


plin ^ (MpX,. - u^x^ 
111 / 


W 




(16) 
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This corresponds to setting 




e 


(17) 


in (d in the linear limit. In the rest frame Uj = 0, ito = 1, we have W = |x'| = r, and then 
the nonvanishing components of Eq. (ITS)) constitute the Coulomb electric field of a point-like 
charge: 


TT’/lin 

^Oi 


x', e 

r dTrr^ 


Expression (I16|) satisfies, in the moving frame, the linear Maxwell equation 


(18) 


Q/J plin 



(19) 


with the current (jB)) of a point charge. This rather evident statement is explicitly demon¬ 
strated in E 

In order to reproduce the standard form of the electric and magnetic field components 
produced by a charge moving with the time-independent speed = 0 with its worldline 
X = passing through the origin x = 0 at zero time-moment xq = 0, written in Ref. [J 
in the case of linear Maxwell electrodynamics as 


Eiin 


e(x-^xo) 
c2 ’ 47ri?*3 


[v X E*'"], 


it is sufficient to note that R* defined in [3| (when v has only the first component, Vi 
as 


( 20 ) 


kiv) 


= {xi-^xof + {xl + xl) ( 21 ) 

is related to the Lorentz scalar W (fT2ll in the following way R*^ = ^1 — (Note that 

R' of Ref. 0 is just our invariant W). With this substitution Eq. (EIH) is just what follows 
for the electric and magnetic components of (ITHll . 

The electric and magnetic fields (uni) or m can be written also as functions of the 
difference Aa;^ = — x^ between the coordinate of the observation point x^ and that of 

the charge x^ as follows 




Ri (Aa;,f) = eijkPf^ = eijk 


e {u^Ax^ — u^Ax^) 

47rIE2 W ’ 

e {ujAxk — UkAxj) 


AttW^ 


W 


with 


= 1 - 


-1 


(Axi-^Axof + ({Ax2f + (Axsf^ ^1 - 


( 22 ) 


(23) 


As long as the trajectory is fixed, neither (1^ . nor (1^ contain the variable explicitly. 

In order to write the fields in the Lienard-Wiechert form it is necessary to exploit the 
light-cone condition (Aa:)^ = 0, i.e. (Aa: 2 )^ -b (Axs)'^ = (Axq)^ — (Axi)^ in (1^ . which tells 
that the charge and the observation point must be separated by a light-like interval for the 
field produced by the charge be nonzero in the observation point. Imposing the light-cone 
condition turns (1231) into 


= 
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and (l2^ becomes Eqs. (63,8), (63,9) of Ref. 0 with v = 0 after the identihcation R = Axq, 

Ri — Axi- 

In this connection, the following remark is in order. In linear electrodynamics the influence 
of a charge propagates with the speed of light in the vacuum c, and hence the observation 
point must be separated from the four-position of the source by a light-like interval, (Ax)'^ = 
0, other space-time points carrying no field produced by the source at this position. This 
is not the case in nonlinear electrodynamics. It is well known already in QED that the 
nonlinearity leads Q to nontrivial dielectric permeability and magnetic susceptibility of the 
vacuum in an external field, thereby to deviation of the speed of propagation from that of 
light. The role of the external field is in our case played by the field produced by the charge 
itself, the propagation speed depending, as a matter of fact, on its intensity. Eor this reason 
it will be not relevant to impose the light-cone condition onto the nonlinear solution of the 
next section for getting an analog of the Lienard-Wiechert potential. 


C. Solution to nonlinear field equations 


From (jH) we can conclude that its solution may be related to given by (fTBl) as 

(I - ad{x))F>^'' = (24) 

since (fT51) satisfies the linear Maxwell equation (fTOl) as demonstrated in |Vl Then, using the 
ansatz da and taking into account that = 1, we get 


This is the cubic equation for g{W^): 

3 25 

g3 , 


(l-a;?(x))g = ^^, 
Six) = 


2 e 


aW^ 47raW^ 
The discriminant of this equation is positive 


= 0 . 


^ (saW^) “^(dTrairO ^ 


because > 0, as argued before. The only real Cardano solution for the case is 


g = 


4TraW^ 

Its asymptotic behavior near the charge 

/ 2e 


\/Q + ^ 

\ 1/3 


4t:(jW^ 


-Vq- 




(25) 


(26) 


(27) 


(28) 


(29) 


will provide convergence for the field mass integral in Subsection ’’Field mass”. In the 
rest frame, when substituted into m with uo = l, Ui=0,W = r this reproduces the 
result of Ref. Q. In the limit of vanishing nonlinearity cr —)> 0 Eq. (^51) becomes g = 
therefore the solution Eq. da turns into expression (1161) for the known field of a moving 


T’lin 


charge in the linear Maxwell electrodynamics, F^^, = F^^^. 

Finally, from da for the point charge moving along the trajectory x = with a 
constant speed along axis 1, expressions for the electric and magnetic fields Ei = and 
Bi = eijkFjk as functions of the observation coordinates are 


E = (1 - !^)-i/2(x-Ia:o)5(IE'), B = [v x E]. 


(30) 
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In terms of the distance from the charge Aa;^ = these are 

/ 2 \ “1/2 

E=(^l-y 5(W^")(Ax-^Aa:o), 

/ 2\ “1/2 / 2\ “1/2 
Si=0, B2 = -2^h-^j 5(W"")Aa:3, B3(Ax) = 2^1 - ^j 5(W"")Ax2. 

(31) 

where is the same as (1251) . and g{W^) is the solution (E51) . This is the nonlinear 
generalization of Eqs. (uni) and (1^^ . 


IV. SOLITONIC REPRESENTATION OF THE POINT-CHARGE FIELD 


A. Energy-momentum tensor 

In this section we set c = 1, and also = 1 in the expression for a. The Noether current 
corresponding to space-time translations 

^ dA^ dLjx) _ 

dx^d{dA>'/dx^) ^ ^ 

calculated using the Lagrange density © with the account of the field equations o is 
written as 

= - [1 - ^ [A>^ (1 - a;?(x)) Fj)] - g^^''L{x). (32) 

This is the energy-momentum tensor of the electromagnetic field. The requirement that it 
should be gauge-invariant makes us omit the potential-dependent second term (that is a full 
derivative) to be left with the expression, symmetrical under the permutation of the indices: 


= - [1 - crS^(a:)] F'^^Fj) - g^''L{x). 

Densities of the energy and of the Pointing vector have the form 


(33) 


r°° = -F°"F° (1 - al?(x)) +l?(x)-A (5(x))^ , 


3 OA pQ , 
^i\ 77.0 , 


pM ^ _p^^p0 ^ 

Contervariant components of the tensor (13311 are 

Tl^ = - [1 - a;?(a;)] F'^^F^a - Si^Lix). 

Its trace is different from zero 

T: = 2ai:S{x)f . 


(34) 

(35) 

(36) 


in the nonlinear case cr 0. 


1. Weak continuity and energy-momentum conservation 

The above construction of the energy-momentum tensor does not provide the fulhllment 
of its continuity, since the current is meant to be supported by outer forces and therefore we 
should admit the energy-momentum non-conservation. Instead of the continuity equation 
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we formally obtain what may be referred to as the ’’partial conservation of the Noether 
current” in the form 


dx>^ 


Ti; = 


(37) 


To derive this relation, equation of motion (jdj) was used in Appendix 2 together with the 
Bianchi identities This relation is not owing to the nonlinearity and 

retains its form in the standard linear theory as well. The right-hand side in (I37p is not 
vanishing. However the continuity of the energy-momentum tensor holds in a weak form 


= ( 38 ) 

To see this, note that the 4-current of a point charge moving with the constant speed is 
parallel to its 4-velocity vector. Then the right-hand side of (I37p disappears when contracted 
with the 4-velocity, u'^Fj^xj^ ^ u'^Fj^xu^ = 0. The property of weak continuity (155)) will be 
sufficient for establishing the conservation of the energy-momentum vector. 

Let us define the latter by the integral over a space-like hyperplane that is orthogonal to 
the four-velocity and crosses the time axis at xq = s\/l — 


P^^ = 


/ 


u'^TH’d (ux 


s) d^x. 


(39) 


Via the Gauss theorem it follows from the vanishing of the 4-divergence (l38l) that is 
independent of s, and thereby of the time of observation xq, because the hyperplane can be 
shifted as a whole along the vector u'^ without affecting the value of the integral, since the 
fields (fTH) decrease at space-time infinity no less fast as in the linear electrodynamics, the 
nonlinearity fades away far from the charge, where its fields are weak. On the contrary, one 
cannot change to a space-like hyperplane inclined differently than in ()39p in the definition 
of the energy-momentum 4-vector due to the lack of the continuity law analogous to (l38)l 
with the unit vector u'^ other than the 4-velocity. 

When the field F^x, on which the energy-momentum tensor Tjf depends, is that of a 
uniformly moving (or resting) point charge, the integral in (1551) usually diverges and hence 
makes no sense. This is not the case in the nonlinear theory under consideration here, as 
we shall see in the next subsection. Therefore, we may treat the energy-momentum vector 
dSD) seriously. 

Bearing in mind that s is a Lorentz scalar (moreover, set equal to zero in what follows), 
and that (1551) is a tensor, the integral ((551) does define a Minkowski vector. It cannot 
help being directed along since this is the only external vector in the integrand of (1551) . 
Hence we write 


where Mf is the field mass. 




(40) 


B. Finite field mass 


It follows from (155)) that 

Mf = P’^Uv = J UvT^''ufj,6 {ux) d'^a;. 
From (1551) and ()T6l) we calculate 

Mf = J [— [1 — a^{x)]Ufj,F^^FxU^ — L] 6 (ux) d‘^x = 


1 — a 


2 ctIVV 


W^g 


6 {ux) d'^x. 


(41) 
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Thanks to the delta-function we may set = —x^ in the integrand, the argument of the 
function g{W^) included. For any function provided the integrals below converge, 

the following chain of relations holds 

/ {ux) d^a; = — / — x\ — X 3 )da;idx 2 dx 3 . 

J Uo J 

Here we have integrated over xq using the delta-function. Performing the change of the 

.y! 11 ^ _ ! 

variable xi ^ ^ = a;ivl — = x'l and omitting the prime afterwards we find that this 

integral is equal to 


J $(-x2)d^a;. 

Applying this result to (l4T|l we get for the field mass the integral over the space 


Mf = 


^9{^r 3 


+ o (x") 


d^x. 


Note that the function g (x^) involved here is just (l28l) taken in the rest frame. 
With the help of Eos. dM)) the second term can be expressed as 


,4 3 




Then we can rewrite expression for the mass as 


eg (x^) 

47 r|x| 


Mf = 


g(x^)^ 3eg{x‘^) 


+ 


d'^x. 


4 ' 4 47r|x 

Therefore, to calculate the mass Mf, we have to calculate two integrals. The first one is 


(42) 


^g{^^Yd^x=\e\i 




^ 2cr(47r)^ ^ 

CXD _ _ 2 

= y y® ^+ y" - dy = 0.885, 


(43) 


and the second one is 


^d3x=|e|i 

47 r|x| ' ' 


2 cr(47r) 


oo 2 

^2 = y y"^ \/rT^-b 1 - + dy = 3.984. 


(44) 


The final result is 


Mf = |e|2 


2CT(47r)^ 


3/2 - 2^1 ) = 2-65|e|5 


2CT(47r)^ 


2 < 


and we can see that this is the same value as the full electrostatic energy of a stationary 
pointlike particle found in Q. 
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V. CONCLUSION 

The object of the study in this paper is a nonlinear electrodynamics. We represented a 
moving pointlike electric charge as a soliton, a particle-like solution for its electromagnetic 
field that makes up a field configuration with finite energy. The property of finiteness of the 
field mass is inherent in many nonlinear models, for example in ones in Refs. [l^ . and in 
the first place, in the famous Born-Infeld model m- This property has been encouraging 
attempts to attribute the experimental value of the electron mass to the ener gy c arried by 
its field, the electromagnetic contribution being sometimes almost exhausting [12| . 

To specialize our consideration, we have chosen the model of Ref. Q. This nonlinear 
model originates from truncation of quantum electrodynamics at the second power of the 
field invariant jj,and it may be thought of as the simplest model among those that result 
in finiteness of the field energy of a point charge, because it admits analytical solution to 
the held equations. Besides, its Lagrangian © does not possess the disadvantage of being 
a singular function of the held. It is this quality that led to presence of a maximum value of 
electric held, responsible, in the end, for the hniteness of the held energy in the Born-Infeld 
[m, and other [T^, [l^, like models. On the opposite, in the model of the held of a 
point charge is still singular near the charge (1^ . but this singularity is suppressed to the 
extent sufficient for the convergence of the held-mass integral (1411) . 

By explicitly solving the nonlinear Maxwell equation ([4]) in an inertial Lorentz reference 
frame with the use of covariant ansatz m that includes the 4-vector of the particle 
speed, we found the electric and magnetic helds of a uniformly moving charge as functions 
of the observation point (l30)) and of the 4-distance to the moving charge (1311) . 

We dehned a gauge-invariant symmetric energy-momentum tensor (1331) for the held of 
a moving charge. Irrespective of a special choice of nonlinear model and already in the 
linear limit cr = 0, the energy-momentum tensor possesses but a weak conservation property 
dSZl), dMl), because the current, corresponding to the charge moving with a constant speed 
is external. We were able to establish, nevertheless, the time-independence of the energy 
and momentum as components of the Minkowskian vector dehned as an integral p9D of 
the energy-momentum tensor over a hyperplane orthogonal to The energy-momentum 
vector dMl) is the same as the mechanical energy-momentum of a mere massive particle 
(soliton) with the mass equal to the hnite rest energy of the electric held. 

The consideration in the present article can be readily extended to include any local 
nonlinear electrodynamics with convergent held energy. This condition is met if the corre¬ 
sponding Lagrangian in place of o grows with the held invariant already as with w > | 
(to be published elsewhere). 
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Appendix 1 

Here we demonstrate that Eq. m satishes linear equation m with the point-like charge 
current ([8]). We omit the superscript ”lin” within this Appendix. 

It is convenient to rewrite m in terms of coordinates of the coordinate system K ', which 
moves along the axis xi with constant velocity v (rest frame of the moving particle). The 
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corresponding Lorentz boost is, 

a^o = 7 (x'o + 1 xi =j 


x'2 = X2, ^3 = 0:3. 


and the inverse transformation is given by ©■ The Lorentz transformation reduces the 
scalar to the form 


= (xuy -x^ = 

= x'l + X2 + 2 : 3 ^ = r'^. 


(45) 


The relations 


dxQ 

dx[ 

dxo 

dx'o 

dx^ 


dx[ 

dxi 


= 7 


= 7 


= 7 


= 7 


d (xo - ^xi) 


dxo 

d {xi - f Xo) 


dxo 

d (xq - f xi) 


dxi 

d (xi - |xo) 


9X1 


= 7. 


V 

^ c ’ 


= -7- 


= 7, 


(46) 


where 1,^, designates derivative calculated at constant x, follow from Q. Then 


d 

dxo 

d 

dxi 


= 7 


A. 

dx\ 


0 x' 

V d 


V d 
c dx[ 


dx[ 


(47) 


Now, with the account of (H5]) - (IT7l) . one can calculate divergency of exploiting the 
covariant form (HU). The zeroth and first components are: 




47r 

ey 

47r 


d^F^, = 


47r 

V ey 

c 47r 


d /-7(xi-fxo) 

9x1 \ 

9 /-x: 

9x' \ r'3 

9 /y (xi - fXo) 
9xo I 


W3 


9 /- 
9x' 


r'3 


The second and third components are 
9'^F^2.3 = 


+ 


9 /-yx2 


9X2 


9 

9x2 


1T3 


-y^a:2 

1^3 


+ 


This is zero, because 


9 

9x0 


u 9 
c 9xi 


r'-3 = 


V 9 
9x 


7+7 


u 9 
c 9x'i 


„'-3 _ 


= 0 , 


9 

9x3 


9 

9x3 


672 : 2.3 


r 1 ^ 


f ^ 'll 

672 :;, 3 

■ 9 

47r 

cal 

0 1 


c 9xi 


47r 

.9x0 


- 72:3 

1^3 


-772:3 

t^3 


u 9 
c 9xi 


9r' 


,/-3 


dx'o 


= 0 . 


(48) 


- x'—const 
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Referring to the linear Maxwell equation in the rest frame 

il; (^) = 4,^3,x'), (49) 

we can write 

= An (ejS^{x'), ^e7(5^(x'), 0, 0^ . 

This coincides with the current 

Appendix 2 

In this Appendix we present a detailed derivation of the partial conservation law Eq. (EZl) 
for a more general Lagrangian than © 

L(x) = -Six) +£, £ = ^iSix)), (50) 

where £ is an arbitrary function of Six). The second pair of nonlinear Maxwell equations 
in place of (|3]) now reads 


d 

dx'' 



pfiv 


and the stress-energy tensor ((55|l becomes 



(51) 


^ _p>.Xpu (^1-^1^) + g^^'^Six) - 


(52) 


To obtain the partial continuity equation (IWll one needs to calculate the derivative 
on solutions of the field equations, i.e. by using ([5T1) and the Bianchi identities (US]). We 
have 


d 

dx>^ 


- 1 - 


dZ 
dSix) J 


1 




where 




Let us calculate the derivative 


dSjx) ^ IpXp dFxp 
dx"' 2 dx'' ■ 


Next we use HID to transform the derivative as 

dFxp 


Now Eq. (IMl) can be written as 

-«^[-sw+2| = |i- 


dFp, 

dF,x 

dx^ 

dxP 

dZ \ 

1 p>'P 

dSix)} 

2 


dx^ ’ 


dFp, dF,x 

dx^ dxP 


(53) 


(54) 


(55) 


(56) 


Taking into account that F^p is antisymmetric, and making the change of the summation 
indices A —>■ /i in the first term, and p —>■ /x in the second, we can see that the both terms in 
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the latter expression coincide and their sum is just the second term in (1551) with the opposite 
sign to cancel its contribution into the full derivative of T^. Gathering together (1551) and 
we finally obtain dSll). 
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